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What  are  the  major  goals  of  the  project? 

List  the  major  goals  of  the  project  as  stated  in  the  approved  application  or  as  approved  by  the  agency.  If 
the  application  lists  milestones/target  dates  for  important  activities  or  phases  of  the  project,  identify  these 
dates  and  show  actual  completion  dates  or  the  percentage  of  completion.  Generally,  the  goals  will  not 
change  from  one  reporting  period  to  the  next.  However,  if  the  awarding  agency  approved  changes  to  the 
goals  during  the  reporting  period,  list  the  revised  goals  and  objectives.  Also  explain  any  significant 
changes  in  approach  or  methods  from  the  agency  approved  appljcatjon  or  plan. _ 


In  this  work  the  grantees  use  a  combination  of  experimental  and  computational  tools  to  design 
and  produce  compositionally  graded  metallic  glass  alloys  with  desirable  properties,  particularly 
for  the  extreme  conditions  of  ballistic  impact.  This  research  is  in  direct  response  to  Topic  Per5- 
G:  “Amorphous  Metals  for  High-Performance  Penetration”  of  the  solicitation  this  project  re¬ 
sponded  to.  This  topic  is  defined  within  DTRA  Basic  Research  Thrust  Area  4:  “Science  to  De¬ 
feat  WMD.”  The  fabrication  of  metallic  glasses  with  graded  composition  pursued  in  the  current 
project  will  enable  the  manufacture  of  structures  with  graded  properties  for  improved  penetration 
and  payload  delivery,  such  as  an  outer  shell  that  is  hard  and  remains  smooth  on  the  outer  sur¬ 
face  but  tough  in  the  interior. 

The  scope  of  the  overall  project  encompasses  three  technical  objectives: 

•  Proof  of  concept  -  demonstrate  that  graded  metallic  glass  coatings  influence  the  trans¬ 
mission  of  high-kinetic  energy  penetrators  through  sand  layers  positively 

•  Show  that  graded  metallic  glass  layers  can  indeed  be  fabricated,  along  with  identification 
of  their  composition-dependent  glass  forming  ability  and  mechanical  properties 

•  Explore  possible  compositions  and  composition  gradients  through  computations 

On  the  technical  side,  these  three  major  milestones  are  achieved  by  the  following  innovations: 

•  Develop  a  finite-element  model  that  allows  to  explore  penetrating  properties  of  projec¬ 
tiles  with  different  surface  layers,  especially  those  with  graded  properties. 

•  Use  computational  modeling  to  establish  novel,  physics-based  predictive  rules  for  glass 
forming  ability  in  systems  of  2  or  3  elements.  Produce  small  batches  of  candidate  alloy 
compositions  for  rapid  model  validation,  using  a  laser  deposition  technique  to  permit  rap¬ 
id  and  efficient  evaluation  of  the  alloy  structure  and  glass  stability. 

•  Using  the  simple  2  or  3  element  alloys  identified  by  the  computation-based  optimization 
process  as  a  starting  point,  produce  compositionally  graded  alloys  of  multicomponent 
(>3  component)  systems  via  laser  deposition.  Use  this  combinatorial  technique  to  refine 
the  alloy  compositions  for  optimal  glass  forming  ability  and  other  properties.  Produce  the 
most  promising  compositions  in  bulk  via  casting  techniques  for  further  characterization. 
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Year  5  of  the  project  should  have  concentrated  on  designing,  fabricating,  and  testing  alloys  with 
desirable  properties,  such  as  high  strength  under  dynamic  loading,  based  on  the  previous  work. 
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What  was  accomplished  under  these  goals? 

For  this  reporting  period  describe:  1)  major  activities;  2)  specific  objectives;  3)  significant  results,  including 
major  findings,  developments,  or  conclusions  (both  positive  and  negative);  and  4)  key  outcomes  or  other 
achievements.  Include  a  discussion  of  stated  goals  not  met.  As  the  project  progresses,  the  emphasis  in 
reporting  in  this  section  should  shift  from  reporting  activities  to  reporting  accomplishments. 


1 .  Proof  of  concept  through  Continuum  Simulations 

1.1  Summary 

Finite  element  (FE)  simulations  were  performed  to  assess  the  hypothesis  that  metallic 
glass  projectiles  with  graded  properties  would  perform  better  as  penetrators  than  homogenous 
projectiles.  In  the  simulation  a  metallic  glass  projectile  impacts  and  travels  through  dry  sand. 
The  kinetic  energy  and  damage  to  the  projectile  are  tracked  as  a  function  of  depth.  A  second  set 
of  simulations  were  performed  where  after  traveling  through  the  sand  to  a  specified  depth  the 
projectile  impacts  and  penetrates  a  reinforced  concrete  slab.  Two  different  grading  schemes 
were  considered  i)  grading  of  properties  from  front  to  rear  of  the  projectile  and  ii)  and  radial 
grading  from  properties  from  inner  core  to  outer  skin.  Projectiles  with  homogenous  properties 
were  also  explored.  As  expected  in  the  homogenous  case  the  penetration  depth  was  positively 
correlated  with  increasing  density,  elastic  stiffness,  and  compressive  strength.  For  the  graded 
projectiles  grading  density  along  the  length  (from  highest  density  at  the  nose  to  lowest  density  at 
the  rear)  was  found  to  produce  ~9%  increase  in  penetration  depth  over  the  homogenous  projec¬ 
tile.  Projectiles  with  graded  stiffness  and  strength  along  the  length  and  projectile  with  radial 
grading  of  properties  showed  significantly  smaller  improvement  in  penetration  depth. 

1.2  Simulation  Overview 

This  section  serves  as  a  brief  overview  of  the  FE  simulations.  A  complete  description  of 
the  material  properties  can  be  found  in  the  appendix. 


All  simulations  were 
performed  using  a  coupled 
Eulerian/Lagrangian  (CEL) 
approach  in  ABAQUS  explicit. 
The  sand  was  modeled  as  a 
Euleraian  fluid  and  the  projec¬ 
tile  and  concrete  slab  with  an 
standard  Lagrangian.  The 
Eulerian  part  of  the  FE  mesh 
in  CEL  analyses  is  represents 
a  stationary  cube  which  the 
Eulerian  material  (sand)  can 
move  into  and  out  of  as  it  in¬ 
teracts  with  the  Lagrangian 
structures.  The  element  may  be  completely  or  partially  filled,  and  is  tracked  as  it  flows  through 
the  mesh  as  a  local  volume  fraction  in  each  grid.  The  geometry  of  the  simulation  is  shown  in 
figure  1 .  The  Eulerian  sand  is  constrained  in  a  “tube”  of  concrete  to  prevent  loss  of  sand  materi¬ 
al  from  the  sides  of  the  simulation  volume  and  to  mimic  the  constraining  pressure  a  large  body 


Figure  1:  Geometry  of  CEL  FE  simulation.  Sand  is  modeled 
as  an  Eulerain  fluid  while  the  concrete  and  projectile  as  La¬ 
grangian  solids. 
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of  sand  (semi-infinite)  would  apply  in  a  small  simulation  volume.  The  projectile  is  given  an  initial 
velocity  of  600  m/s  in  air  before  impacting  the  sand. 

Zr55Al3oNi5Cu3o  was  chosen  as  the  model  metallic  glass  for  the  projectile  as  experimental 
values  for  the  required  dynamic  material  properties  were  available  in  the  literature  (Mashimo  et 
al  2006).  The  material  was  assumed  to  be  elastically  isotropic  with  a  Youngs  Modulus  of  86GPa 
and  a  shear  modulus  of  31  GPa.  The  yield  surface  was  described  by  a  Drucker-Prager  model.  A 
Mie-Gruneisen  equation  of  state  (EOS)  calibrated  to  shock  Hugoniot  compression  data  was 
used  to  model  the  volumetric  strains  under  dynamic  conditions  and  high  pressures.  Hardening 
was  handled  using  a  Johnson-Cook  plasticity  model  coupled  with  a  progressive  damage  model. 
Details  can  be  found  in  the  appendix. 

The  dynamic  properties  of  sand  were  taken  from  the  DOE  report  “Shock  Response  of 
Dry  Sand”  (Brown  et  al  2007).  An  EOS  model  with  compaction  was  used  to  model  the  defor¬ 
mation  of  sand  under  dynamic  loading.  An  initial  porosity  of  40%  was  assumed  based  on  Brown 
et  al.  As  a  bulk  volume  the  sand  behaves  as  an  elastic  solid  with  a  shear  modulus  of  50  MPa. 
The  sand  begins  to  flow  plastically  at  a  pressure  of  0.9  MPa  and  begin  to  compact.  At  1.1  GPa 
the  sand  is  completely  compacted  and  all  pores  are  crushed. 

The  concrete  properties  were  taken  from  the  numerous  concrete  tutorials  in  the 
ABAQUS  Example  Manual.  These  properties  were  calibrated  to  samples  taken  from  the  Koyna 
dam  and  serve  as  an  archetype  material  for  simulation  of  concrete  structures  under  seismic 
loading.  Damage  in  the  concrete  was  modeling  using  the  Concrete  Damage  Plasticity  model. 
Note  that  this  is  a  simplistic  damage  model  and  the  results  presented  below  of  projec¬ 
tile/concrete  slab  interaction  should  not  be  taken  quantitatively.  Instead  the  simulation  was  in¬ 
tended  to  show  trends  in  impact  energies  as  a  function  of  slab  depth.  Figure  2  shows  a  se¬ 
quence  of  still  images  showing  the  impact  and  penetration  of  the  concrete  slab  by  a  metallic 
glass  projectile. 


Figure  2:  Impact  and  penetration  of  a  concrete  slab  by  a  metallic  glass  projectile.  The  sand  has 
been  removed  for  clarity. 


1.3  Simulation  Results 

1.3.1.  Homogenous  Metallic  Glass  Projectile.  In  order  to  isolate  the  effect  of  individual  prop¬ 
erties  on  penetration  depth  specific  properties  were  perturbed  from  those  measured  for 
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Zr55AI3oNi5Cu3o  metallic  glass.  Properties  were  perturbed  independently  keeping  the  remainder 
of  the  properties  unchanged.  The  penetration  depth  of  the  projectile  in  sand  was  strongly  corre¬ 
lated  with  increasing  density,  elastic  stiffness,  and  compressive  strength  of  the  glass  with  densi¬ 
ty  having  the  largest  effect.  Figure  3  shows  the  main  results. 

1.3.2.  Graded  Metallic  Glass  Projectile.  Two  different  grading  schemes  were  considered  i) 
grading  of  properties  from  front  to  rear  of  the  projectile  and  ii)  and  radial  grading  from  properties 
from  inner  core  to  outer  skin.  The  specific  geometries  are  shown  in  Figure  4. 

Table  1  Properties  for  the  projectiles  with  along  the  length. 

Section _ Density _ Stiffness _ Strength _ 

1  7200  kg/m3  (+3%)  33.0  GPa  (+5%)  1.85  GPa  (+5%) 

2  6790  kg/m3  31.3  GPa  1.77  GPa 

3  6500  kg/m3  (-3%)  39.0  GPa  (+5%)  1.65  GPa  (-7%) 

Table  1  details  the  properties  used  to  explore  the  effect  of  gradients  along  the  length. 
For  the  graded  projectiles  grading  density  along  the  length  (from  highest  density  at  the  nose  to 
lowest  density  at  the  rear)  was  found  to  produce  ~9%  increase  in  penetration  depth.  Reversing 
the  gradient  direction  (lowest  in  nose)  decreased  penetration  depth  by  5%.  Grading  the  stiffness 
along  the  length  (high  to  low)  showed  a  modest  improvement  of  2%.  When  the  stiffness  gradient 
was  reversed  the  penetration  depth  was  virtually  identical  to  the  homogeneous  projectile. 
Strength  was  an  interesting  case  in  that  both  positively  and  negatively  graded  structures  per¬ 
formed  better  than  the  homogenous  glass.  With  an  8%  increase  in  penetration  depth  for  the 
high  to  low  graded  structure  and  a  3%  improvement  for  the  low  to  high.  Figure  5  and  Table  2 
show  the  key  results. 

For  the  radially  graded  structures  no  significant  improvement  or  debit  of  performance 
was  observed.  Results  are  shown  in  the  Appendix. 
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Figure  3:  Trends  in  penetration  depth  for  metallic  glass  projectiles 
with  varying  material  properties. 


Table  2:  Penetration  depth  and  improvement  for  graded  metallic  glass  projectiles. 

Original  Density+  Density-  Stiffness+  Stiffness-  Strength+  Strength- 


Penetration 

(m) 

Improvement 

(%) 


38.93 


42.73 


8.9% 


37.06  39.83 

(4.8%)  2.3% 


39.19  42.11 

0.6%  8.2% 


40.21 


3.3% 
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Material:  Metallic  Glass  (Zr55Al1(}N i5Cu3(j) 
(Lagrangian  Part) 

1.  Homogenous  Section 
Radius:  0.2 
Length:  1 

Round:  0.2 

2.  Graded  properties  along  length 

The  projectile  is  divided  in  to  three  parts  along 
length.  Material  properties  are  varied  in  each 
layer. 

3.  Graded  properties  along  thickness 

The  projectile  is  extruded  by  an  inner  offset  part 
and  compensate  with  a  same  size  part  with 
different  properties. 


Figure  4:  Geometry  of  the  graded  projectiles. 


7 


Figure  5:  Velocity  vs  penetration  depth  for  graded  projectiles  as  compared  to  a  homogenous 
metallic  glass  projectile.  +  indicates  gradients  as  given  by  table  1  and  -  indicates  the  opposite 
direction  (3  high  and  1  low). 


2.  Fabrication  and  Testing  of  Graded  Metallic  Glass  Layers 

2.1  Summary:  Laser  deposition  as  a  combinatorial  BMG  design  tool 

While  unresolved  issues  prevent  the  realization  of  fully  amorphous  3-dimensional  BMG 
parts  prepared  by  LD,1'4  single  layer  deposits  nearly  void  of  detectable  crystalline  features  have 
been  successfully  demonstrated  by  co-PI  Flores.5  In  addition,  several  groups  have  reported 
building  compositionally  graded  crystalline  specimens  with  LD.6"12  Combining  the  two  approach¬ 
es,  this  project  had  as  a  focal  goal  to  establish  LD  as  a  promising  combinatorial  tool  for  develop¬ 
ing  novel  BMGs  with  desirable  properties,  where  both  composition-dependence  and  glass¬ 
forming  ability  can  be  studied  by  varying  powder-feed  rates  and  laser  power  in  single  or  few-run 
depositions.  As  a  proof-of-concept  and  method-development  study,  we  examined  the  glass 
forming  Cu^Zr*  binary  system  over  a  wide  composition  range  and  confirmed  the  previously  re¬ 
ported  main  glass  forming  compositions  of  Cu65Zr35  and  Cu5oZr5o. 

For  the  ternary  Cu-Zr-Ti  system,  we  could  equally  show  that  our  fabrication  scheme  al¬ 
lows  to  systematically  probe  large  regions  of  ternary  composition  space.  We  identified 
Cu51.7Zr36.7Tiii.6  as  the  best  glass  former  in  a  fraction  of  the  experimental  time  of  previous  study, 
which  identified  similar  compositions  of  Cu57.6Zr32.4Tho, 13  and  Cu60Zr3oTi10.14 

Finally,  we  have  demonstrated  that  nanoindentation  along  with  appropriate  analysis  can 
be  used  on  the  graded  LD  samples  to  map  out  mechanical  properties  with  high  accuracy,  thus 
completing  a  high-throughput  capability  to  identify  glass  forming  alloys  and  their  mechanical 
properties  from  a  large  composition  space. 

2.2  Continuously  graded  Cu-Zr 

As  a  prelude  to  developing  novel  multicomponent  (3  or  more  elements)  BMG  alloys,  the 
first  experimental  task  was  to  validate  the  combinatorial  methodology  by  fabricating  graded  bi¬ 
nary  specimens  containing  known  BMG-forming  compositions.  Among  binary  alloys,  the  exten¬ 
sively  studied  Cu-Zr  system  was  the  obvious  choice,  due  to  the  existence  of  several  composi¬ 
tions  with  reported  casting  thicknesses  of  1  mm  or  larger.15,16 

All  previously  reported  efforts  to  fabricate  compositionally  graded  specimens  were  unsuccessful 
and  varied  the  composition  of  the  build  from  layer  to  layer  by  changing  the  relative  delivery  rates 
of  two  different  powders.  In  order  to  overcome  these  problems,  we  have  developed  a  two- 
dimensional  approach  to  fabricate  single-layer  specimens  that  feature  a  compositional  gradient 
across  the  surface  area  of  the  deposit. 

2.2.1  Fabricating  continuously  graded  Cu-Zr  specimens  by  LD 

All  LD  experiments  were  carried  out  with  an  Optomec  MR-7  Laser  Engineered  Net 
Shaping  (LENS™)  system.  Copper  and  zirconium  feedstock  powders  with  minimum  purities  of 
99.9  and  99.2  at.  %,  respectively,  were  deposited  on  grade  702  Zr  substrates  with  approximate 
dimensions  of  3.5  cm  x  4.0  cm  x  0.5  cm.  Prior  to  deposition,  the  substrates  were  surface-ground 
with  320  grit  SiC  paper  and  cleaned  with  methanol  to  remove  any  surface  contaminants.  To  fab¬ 
ricate  the  graded  Cu-Zr  specimens,  a  square  deposition  geometry  was  chosen  with  dimensions 
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of  25.4  x  25.4  mm.  Cu  and  Zr  powders  were  simultaneously  deposited  with  a  250  W  laser  trav¬ 
eling  at  8.5  mm/s.  To  achieve  compositional  gradation,  the  delivery  rate  of  both  powders  was 
varied  in  a  stepwise  manner  from  hatch  to  hatch  along  the  raster  path.  The  Cu  delivery  rate  was 
gradually  increased  from  3.51  to  4.88  g/min  while  the  Zr  delivery  rate  was  decreased  from  3.06 
to  2.36  g/min.  After  the  initial  deposit,  a  200  W  laser  with  a  travel  speed  of  12.7  mm/s  was  used 
to  create  9  re-melted  lines  spaced  2.54  mm  apart,  in  the  direction  of  the  composition  gradient. 
The  re-melting  step  facilitated  adequate  localized  melting  and  mixing  of  the  copper  and  zirconi¬ 
um  powders.  The  alloyed  lines  were  then  re-melted  a  final  time  with  different  laser  powers  rang¬ 
ing  from  100  to  180  W  at  a  higher  travel  speed  of  16.9  mm/s. 

Figure  6  shows  the  compositional  profile  of  the  line  re-melted  by  a  100  W  laser,  meas¬ 
ured  with  energy  dispersive  x-ray  spectroscopy  (EDS).  The  measurements  were  performed  af¬ 
ter  polishing  the  entire  specimen  to  a  mirror  finish  with  colloidal  silica.  The  composition  profile 
was  continuously  graded  from  61-76  at.  %  Cu  and  is  well  described  by  a  linear  trend  line.  The 
predicted  compositional  profile  based  on  the  powder  delivery  rates  was  in  excellent  agreement 
with  the  measured  compositional  trend,  confirming  the  ability  of  the  experimental  methodology 
to  accurately  control  the  compositions  of  graded  specimens  by  choosing  appropriate  powder 
delivery  rates. 

2.2. 1 1dentifying,  yjtrifiecl  regions  in  the  Citrich 
graded_  CitZr  specimens 

In  addition  to  the  capability  to  rapidly 
synthesize  compositional  libraries,  a  complete 
combinatorial  methodology  for  developing  new 
metallic  glass  alloys  must  include  a  means  for 
high-throughput  screening  of  amorphous  mate¬ 
rial  within  the  library.  X-ray  and  electron  diffrac¬ 
tion  are  the  most  dependable  techniques  for 
verifying  amorphous  structure.  However,  con¬ 
ventional  x-ray  diffraction  techniques  are  limited 
to  interrogating  relatively  large  (millimeter-scale) 
areas,  while  micro-XRD  and  electron  diffraction 
often  requires  several  hours  to  prepare  and  ana¬ 
lyze  a  single  composition.  Therefore,  they  are 
not  ideal  for  the  rapid  screening  of  metallic 
glass-forming  regions  in  large  compositional  li¬ 
braries. 

As  an  alternative,  Ding  et  al .21  reported  a 
high-throughput  approach  by  distinguishing 
amorphous  from  crystalline  material  by  inspec¬ 
tion  of  surface  topography.  They  showed  that  for 
unconstrained  solidification  as  we  have  here,  a 
smooth,  featureless  topography  is  indicative  of 
glass  formation,  whereas  crystalline  features  form  visibly  distinct  asperities. 

Preliminary  identification  of  potential  glass-forming  compositions  in  the  compositionally- 
graded  Cu-Zr  specimens  was  accomplished  using  a  Nikon  light  microscope  equipped  with  dif¬ 
ferential  interference  contrast  (DIC)  for  enhanced  imaging  of  the  specimen’s  three-dimensional 
surface  topography.  In  addition,  DIC  microscopy  is  capable  of  detecting  surface  asperities  with 


Position  (mm) 

Figure  6.  Compositional  profile  of  a  Cu-Zr 
deposit  line  re-melted  using  a  laser  power  of 
100  W  and  a  travel  speed  of  16.9  mm/s. 
The  measured  profile  is  compared  with  the 
ideal  profile  predicted  by  the  powder  feed 
rates,  assuming  perfect  deposition  efficien¬ 
cy  of  both  the  Cu  and  Zr  powders. 
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diameters  as  small  as  40 
nm,17  making  it  a  suitable 
high-throughput  characteriza¬ 
tion  tool  for  identifying  amor¬ 
phous  material  in  the  graded 
specimens. 

Figure  7(a)  shows  a 
DIC  image  of  a  possible 
amorphous  segment  within 
the  line  re-melted  by  the  100 
W  laser.  The  surface  of  the 
segment  was  topographically 
smooth  with  very  few  identifi¬ 
able  surface  features.  In  con¬ 
trast  to  the  smooth  segment, 
the  rough  and  uneven  topog¬ 
raphy  of  the  material  to  the 
immediate  left  of  the  segment 
clearly  indicated  the  presence 
of  crystals  in  this  region  (Fig¬ 
ure  7(c)).  To  corroborate  the 
DIC  results,  SEM  micro¬ 
graphs  of  the  polished  speci¬ 
men  obtained  from  the  middle  of  the  segment  appeared  featureless  (Figure  7(b)),  as  should  be 
expected  from  the  homogeneous  microstructure  of  metallic  glass.  Furthermore,  an  electron  mi¬ 
crograph  obtained  near  the  end  of  the  smooth  segment  displayed  crystallites  of  an  unidentified 
phase  embedded  within  a  featureless  matrix  (Figure  7(d)),  suggesting  a  semi-amorphous  transi¬ 
tion  region  between  the  glass  and  crystalline  regions.  EDS  measurements  across  this  transition 
region  yielded  a  compositional  range  of  less  than  1  at.  %.  The  narrowness  of  this  transition  re¬ 
gion  may  reflect  the  acute  sensitivity  of  GFA  to  compositional  deviations  in  the  binary  Cu-Zr  sys¬ 
tem.16 


Figure  7.  (a)  DIC  image  of  a  topographically  smooth  segment 
of  the  compositionally  graded  Cu-Zr  deposit  and  (b)  the  corre¬ 
sponding  secondary  electron  SEM  image  of  the  microstruc- 
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ture,  from  a  line  re-melted  using  a  laser  power  of  100  W  and  Electron  diffraction  in 

travel  speed  of  16.9  mm/s.  The  Cu  content  increases  from  left  the  transmission  electron  mi- 
to  right,  (c)  DIC  image  of  crystalline  material  adjacent  to  the  croscope  (TEM)  was  per- 
smooth  segment,  (d)  The  SEM  micrograph  recorded  from  the  formed  to  verify  the  amor- 
semi-amorphous  region  shows  crystallites  embedded  within  a  phous  structure  of  the  materi- 
featureless  matrix.  al  that  was  identified  by  DIC. 

The  TEM  sample  was  prepared  by  polishing  and  ion  milling  a  strip  of  the  topographically  smooth 
region  from  the  line  re-melted  by  the  100  W  laser.  The  resulting  diffraction  pattern  featured  con¬ 
centric  halos  that  are  characteristic  of  amorphous  materials,  validating  the  reliability  of  DIC  im¬ 
aging  for  identifying  metallic  glass  regions  in  the  graded  specimens. 

An  additional  compositionally-graded  specimen,  also  prepared  using  a  final  re-melt  pow¬ 
er  of  100  W,  was  used  to  investigate  the  glass  transition  and  crystallization  behavior  of  the 
amorphous  region  with  differential  scanning  calorimetry  (DSC).  The  region  was  mechanically 
separated  from  the  substrate  and  heated  at  a  rate  of  20  K/s.  Figure  8  displays  the  heating  trace 
for  the  strip  of  the  featureless  material.  Although  the  glass  transition  was  weak,  the  curve  dis¬ 
played  a  prominent  exothermic  crystallization  peak,  suggesting  that  the  topographically  feature¬ 
less  material  contained  a  significant  volume  fraction  of  metallic  glass.  Tg  was  755  K,  and  the  on¬ 
set  temperature  of  the  crystallization  peak  (Tx)  was  780  K,  consistent  with  the  Tg  and  Tx  of  the 
best  reported  glass  former  in  the  Cu-Zr  system,  Cu64.5Zr35i5.18 


2.2.2  Identifying  the  composition  with  the  highest  GFA 


An  important  parameter  to  consider  during  laser  re-melting  is  the  laser  line  energy  A,  a 
parameter  quantifying  the  instantaneous  heat  input  per  unit  area  by  the  laser  to  the  target  sub¬ 
strate  material.  It  is  defined  mathematically  as  A  =  PI(vD)  where  P  is  the  (absorbed)  laser  power, 
v  is  the  laser  travel  speed,  and  D  is  the  diameter  of  the  laser  beam  on  the  surface  of  the  sub¬ 
strate,  assumed  to  be  independent  of  laser  power  and  set  to  1  mm  in  our  study.  In  laser  deposi¬ 
tion  experiments  with  Zr58.5Cui5.6Ni12.8AI1o.3Nb2.8  powder,  an  excellent  glass  former,  Sun  and  Flo¬ 
res  observed  that  for  a  fixed  laser  power,  reducing  A  below  a  threshold  value  by  increasing  the 


travel  speed  nearly  eliminated  crystalliza¬ 
tion  in  the  heat  affected  zone  (HAZ).5,19  The 
same  effect  of  laser  travel  speed  was  no¬ 
ticed  in  laser  surface  melting  experiments 
of  Cu6oZr3oTi10  substrates.20  Additionally, 
they  demonstrated  that  for  a  fixed  travel 
speed,  decreasing  the  laser  power  below  a 
threshold  value  was  necessary  for  vitrifica¬ 
tion  of  the  substrate’s  surface.  These  ob¬ 
servations  of  maximum  laser  power  and 
minimum  travel  speed  as  criteria  for  vitrifi¬ 
cation  suggest  that  a  maximum  value  of  A 
exists  for  a  given  alloy  composition  ( Amax ), 
above  which  the  thermal  burden  becomes 
too  large  to  bypass  crystallization  during 
the  cooling  of  the  melt  pool. 
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Figure  8.  DSC  heating  trace  of  the  topograph¬ 
ically  smooth  segment  mechanically  separated 
from  the  line  re-melted  using  a  laser  power  of 
100  W  and  travel  speed  of  16.9  mm/s.  DSC  daW 
was  obtained  at  a  heating  rate  of  20  K/min. 


To  investigate  the  effect  of  A  on  our  graded  specimen,  nine  evenly-spaced  re-melted 
lines  were  made  in  the  direction  of  the  composition  gradient,  with  laser  power  increasing  from 
100  to  180  W  in  steps  of  10  W.  Figure  9  shows  a  DIC  image  of  the  complete  specimen.  The  ob¬ 
served  amorphous  segment  in  each  re-melted  line  is  enclosed  by  dashed  lines.  The  width  of  the 
segments  deceases  with  increasing  laser  power  or,  equivalently,  increasing  A.  Although  the 
endpoints  of  the  amorphous  segments  could  not  be  sharply  defined  due  to  the  presence  of  the 
semi-amorphous  transition  regions,  they  were  approximated  from  the  DIC  image  as  the  posi¬ 
tions  where  the  topographically  featureless  matrix  was  no  longer  clearly  visible. 

The  chemical  composition  at  the  two  endpoints  of  each  amorphous  segment  was  meas¬ 
ured  by  EDS,  and  the  results  are  plotted  versus  A  in  Figure  10.  To  achieve  statistical  precision, 
each  point  on  the  plot  represents  the  average  of  five  measurements  at  the  same  position  in  the 
specimen.  The  upper  and  lower  bound  data  are  each  well-fitted  by  linear  trend  lines.  From  the 
intersection  of  these  two  lines  it  may  be  deduced  that  the  optimum  glass  forming  composition  in 
the  graded  specimen  was  Cu64.7Zr35.3,  in  excellent  agreement  with  the  composition  that  has 
been  reported  in  literature,  Cu64.5Zr35.5. 16,18  We  further  deduce  that  the  maximum  laser  line  ener¬ 
gy  for  this  optimal  composition,  above  which  a  glass  is  not  formed,  is  10.9  J/mm2.  While  this 
agreement  is  remarkable,  it  should  be  noted  that  our  predicted  optimal  composition  is  subject  to 
inherent  errors  associated  with  the  accuracy  of  the  EDS  measurement  technique,  as  well  as 
“operator  error”  associated  with  identifying  the  transition  from  amorphous  (featureless)  to  crys¬ 
talline  at  each  endpoint.  We  estimate  that  each  of  these  errors  contribute  ±  1  at.  %  or  less  to 
each  measured  endpoint  and  to  the  predicted  optimal  composition.  If  additional  precision  is  re¬ 
quired,  a  more  gradual  compositional  gradient  may  make  the  amorphous/crystalline  transition 
more  easily  discernable,  reducing  the  contribution  of  “operator  error”. 

The  linear  narrowing  of  the  range  of  glassy  compositions  with  increasing  laser  line  ener¬ 
gy  towards  the  optimum  glass  former  implies  a  direct  relationship  between  GFA  and  Amax.  High 
GFA  is  reflected  in  an  alloy’s  ability  to  carry  a  larger  thermal  burden  without  crystallizing. 


Figure  9.  DIC  image  of  a  graded  Cu-Zr  deposit,  showing 
several  lines  re-melted  at  different  laser  powers.  The  laser 
travel  speed  was  maintained  at  16.9  mm/s.  The  dashi§ 
boxes  indicate  topographically  smooth  segments  of  each 
line,  suggesting  that  the  material  is  amorphous. 


This  ability  to  isolate  and 
identify  the  best  glass  formers  within 
a  composition  library  is  a  marked 
advantage  of  laser  deposition  over 
other  combinatorial  techniques  such 
as  sputtering.  While  smooth  compo¬ 
sition  profiles  could  be  obtained  by 
magnetron  sputtering,  the  excep¬ 
tionally  high  quench  rate  of  sputter¬ 
ing,  and  the  resulting  ease  with 
which  an  amorphous  structure  is 
formed  over  a  broad  composition 
range,  makes  it  impossible  to  direct¬ 
ly  differentiate  one  composition  from 
another  with  respect  to  GFA  without 
relying  on  “predictive”  parameters 
such  as  T^.15,21 

2.2.3  Identifying  additional  Cu-Zr 
glass  formers 

In  addition  to  Cu64.5Zr35.5, 
several  other  compositions  in  the 
Cu-Zr  system  have  been  widely  rec¬ 
ognized  as  bulk  glass  formers  with 
critical  casting  thicknesses  exceed¬ 
ing  1  mm.22,23  Li  and  co-workers 
casted  a  series  of  fifteen  Cu-Zr  wedges  ranging  46-70  at.  %  Cu  and  identified  three  “GFA 
peaks”  corresponding  to  Cu64Zr36,  Cu56Zr44,  and  Cu5oZr5o.15  Although  our  combinatorial  method¬ 
ology  demonstrated  the  accurate  identification  of  the  optimum  glass  former  in  the  binary  system, 
complete  validation  of  the  method  would  necessitate  the  ability  to  identify  all  other  notable  glass 
formers  as  well.  To  thoroughly  investigate  GFA  in  the  Cu-Zr  system,  we  fabricated  continuously 
graded  specimens  with  an  expanded  composition  range. 

The  new  specimen  was  synthesized  according  the  experimental  procedure  that  was 
previously  described  for  the  original  graded  specimen.  However,  during  deposition,  the  range  of 
powder  delivery  rate  for  both  feedstock  powders  was  expanded  to  widen  the  composition  range 
of  the  finished  specimen:  1.78  -  5.24  g/min  and  3.91  -  2.29  g/min  for  the  Cu  and  Zr  powders  re¬ 
spectively.  To  adjust  for  the  lower  GFA  of  the  other  glass  formers,  the  laser  powers  for  the  final 
re-melt  was  reduced  to  range  from  50  to  90  W,  with  increments  of  5  W  for  each  re-melted  line. 
Another  identical  specimen  was  also  fabricated,  albeit  with  laser  powers  ranging  from  50  to  130 
W,  in  steps  of  10  W.  Figure  1 1  provides  a  DIC  image  of  one  of  the  completed  specimens. 

Using  the  described  method  for  deducing  the  best  glass  formers,  an  additional  GFA 
peak  located  at  Cu50.2Zr49.8  was  identified  in  the  extended  range  specimens.  Figure  12(a)  is  a 
plot  of  the  combined  data  for  the  original  specimen  as  well  as  both  extended  range  specimens. 
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Figure  10.  Composition  of  the  graded  Cu-Zr  deposit 
at  the  endpoints  of  the  glassy  regions  observed  via 
DIC  (Figure  4)  and  SEM,  as  measured  using  EDS.  As 
the  laser  line  energy  increases,  the  compositional 
width  of  the  glassy  segment  shrinks  until  only  the  best 
glass-forming  composition  remains. 
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Figure  11.  DIC  image  of  a  graded  Cu-Zr  specimen  with  an  extended  composition  range.  The 
re-melted  lines  in  image  correspond  to  nine  different  laser  powers,  ranging  50  -  90  W,  in  steps 
of  5  W  and  decreasing  from  left  to  right.  The  boxed-in  areas  are  topographically  smooth  seg¬ 
ments  observed  in  the  re-melted  lines.  These  segments  correspond  to  two  GFA  peaks  in  the 
Cu-Zr  binary  system. 
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Figure  12.  (a,  left)  Composition  of  the  graded  Cu-Zr  deposit  at  the  endpoints  of  the  glassy  re¬ 
gions  observed  via  DIC  (Figure  8)  and  SEM,  as  measured  using  EDS.  An  additional  GFA  peak 
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was  observed  at  Cu50.2Zr48.8.  The  overall  results,  agree  well  with  wedge  casting  experiments 
performed  by  Li  et  al .15  (b,  right)  Plot  of  critical  casting  thickness  of  Cu-Zr  binary  alloys.  Image 
is  taken  from  Ref.  1 5. 

For  comparison,  Figure  12(b)  shows  the  GFA  plot  from  Li  et  al.’ s  wedge  casting  experi¬ 
ments.  The  extrapolated  peaks  line  up  nearly  perfectly  with  two  of  the  three  peaks  reported  by 
the  authors  but  the  peak  at  Cu56Zr44  was  not  observed  at  all  in  the  graded  specimens. 

A  factor  inherent  to  the  laser  experimental  procedure  could  account  for  the  discrepancy. 
According  to  the  sequence  of  fabrication  steps,  vitrification  occurs  during  the  final  re-melting 
step,  when  both  a  faster  travel  speed  and  reduced  laser  powers  are  applied  to  increase  the 
subsequent  cooling  rate  within  the  melt  pool.  Prior  to  this  final  step,  the  deposit  is  fully  crystal¬ 
line,  and  depending  on  the  kinetics  of  melting,  certain  crystalline  phase(s)  may  not  fully  melt  at 
low  laser  powers,  leaving  behind  residual  nuclei  that  could  promote  subsequent  heterogeneous 
crystallization  events. 

TEM  samples  from  the  center  and  edges  of  the  glass  forming  segments  in  the  extended 
range  specimens  are  currently  being  prepared  via  the  focused  ion  beam  for  selected  area  dif¬ 
fraction.  The  structural  analysis  will  confirm  the  glassiness  of  the  smooth  regions  and  shed  light 
regarding  the  role  and  effect  of  competing  crystalline  phases  on  GFA. 

2.3  Continuously-Graded  Cu-Zr-Ti  Compositional  Libraries 

Following  the  successful  work  exploring  glass  formation  in  the  binary  Cu-Zr  system,  we 
extended  the  methodology  to  investigate  the  ternary  Cu-Zr-Ti  system.  It  is  known  that  minor 
alloying  of  Cu-Zr  alloys  with  Ti  and/or  Al  could  vastly  improve  the  GFA  and  enhance  mechanical 
properties45-49.  Although  speculation  exists  as  to  the  physical  reasons  minor  alloy  additions  of 
certain  elements  improve  GFA  and  mechanical  properties  of  BMGs,  the  proposed  reasons  are 
mostly  qualitative  and  do  not  provide  a  reliable  means  of  predicting  the  specific  effect  of  a  minor 
alloying  element  in  an  arbitrary  MG  alloy. 

2.3. 1  Fabrication  of  continuously  graded  Cu-Zr-Ti  specimens  by  LD 

An  alternative  fabrication  scheme  was  used  for  the  compositionally  graded  Cu-Zr-Ti  ter¬ 
nary  specimens.  Instead  of  co-depositing  the  Cu  and  Zr  powders  with  varying  powder  delivery 
rates,  three  graded  layers  of  Zr,  Cu,  and  Ti  were  successively  deposited  one  on  top  of  the  other 
in  a  25.4  x  25.4  mm  square  geometry.  The  underlying  substrate  was  the  same  material  as  be¬ 
fore,  grade  702  zirconium  ground  with  320  grit  SiC  paper  and  cleaned  with  methanol  to  remove 
surface  contaminants.  The  minimum  purities  of  the  powders  were  99.9,  99.2,  and  99.2  at.  %  for 
Cu,  Zr,  and  Ti  respectively. 

To  achieve  the  graded  effect  for  a  single  layer  deposit,  the  powder  delivery  rate  was 
fixed  at  a  constant  value  and  the  laser  travel  speed  was  varied  stepwise  from  hatch  to  hatch. 
For  the  initial  layer,  Zr  powder  was  deposited  with  a  250  W  laser,  and  the  travel  speed  was 
ramped  up  from  8.5  to  29.6  mm/s.  Next,  Cu  powder  was  deposited  with  a  250  W  laser  on  top  of 
the  first  layer,  following  the  same  tool  path.  However,  the  travel  speed  was  ramped  down  from 
29.6  to  8.5  mm/s,  resulting  in  the  specimen  being  Zr-rich  on  one  end  and  Cu-rich  on  the  other. 
The  Cu  and  Zr  powder  delivery  rates  were  fixed  at  4.20  and  2.83  g/min  respectively. 

Following  the  two  deposits,  nine  evenly-spaced,  re-melted  lines  were  drawn  across  the 
specimen  with  a  200  W  laser  traveling  at  12.2  mm/s.  The  re-melting  step  was  necessary  to  fa¬ 
cilitate  adequate  mixing  of  the  two  components.  However,  whereas  the  re-melted  lines  in  the 
binary  Cu-Zr  specimen  were  drawn  parallel  to  the  direction  of  the  composition  gradient,  the  re- 
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melted  lines  in  the  case  of  the  ternary  specimen  were  created  perpendicular  to  the  Cu-Zr  com¬ 
position  gradient.  The  aim  for  the  change  in  deposition  direction  was  for  each  of  the  re-melted 
lines  to  hold  a  constant  Cu:Zr  molar  ratio. 

The  final  fabrication  step  was  to  integrate  a  Ti  concentration  gradient  into  each  of  the  re¬ 
melted  lines.  This  was  accomplished  by  depositing  a  layer  of  Ti  in  the  same  manner  as  with  the 
previous  two  layers,  but  using  a  laser  tool  path  that  has  been  rotated  by  ninety  degrees.  The 
powder  delivery  rate  was  fixed  at  0.85  g/min;  the  travel  speed  was  incrementally  changed  from 
8.5  to  29.6  mm/s;  and  the  laser  power  was  set  to  150  W.  After  completion  of  the  third  layer,  the 
nine  re-melted  lines  were  re-melted  again  with  a  200  W  laser  traveling  at  12.7  mm/s  to  thor¬ 
oughly  mix  the  Ti  into  the  alloyed  lines.  A  final  re-melt  of  the  9  lines  with  a  150  W  laser  traveling 
at  a  faster  speed  of  16.9  mm/s  was  performed  to  encourage  vitrification  of  bulk  glass  formers 
within  the  graded  specimen. 


Figure  13.  (a)  Photograph  of  the  complete  ternary  specimen,  (b)  Profile  of  the  XCu/(Xcu+XZr) 
(X  =  atomic  fraction)  molar  ratio  in  the  compositionally  graded  Cu-Zr-Ti  ternary  specimen,  (c)  Ti 
content  in  the  same  ternary  specimen.  The  smooth  contour  fits  demonstrate  the  ability  of  the 
fabrication  scheme  to  systematically  probe  large  regions  of  ternary  composition  space. 
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Figure  14.  DIC  image  of  the  compositionally  graded  Cu-Zr-Ti  specimen.  Three  topographically 
distinct  regions  are  marked,  with  their  corresponding  microstructures  shown  to  the  right  of  the 
large  image.  Regions  1 ,  2,  and  3  are  believed  to  be  amorphous,  semi-amorphous,  and  crystal¬ 
line,  respectively. 

Composition  measurements  were  recorded  by  EDS  across  the  specimen  to  assess  the 
continuity  and  smoothness  of  the  chemical  gradation.  Figures  13b  and  13c  show  plots  of  the 
Xcu/(Xcu+Xzr)  ( X  =  atomic  fraction)  molar  ratio  and  Ti  content  as  a  function  of  position.  Both  sets 
of  data  were  well-fitted  with  smooth  contours,  with  the  direction  of  the  surface  slope  indicating 
that  the  desired  two-dimensional  gradation  had  been  successfully  fabricated  according  to  the 
intended  design.The  DIC  image  of  the  unpolished  specimen  is  shown  in  Figure  14.  Again, 
topographically  smooth  regions  in  the  re-melted  lines  were  easily  distinguished  from  the  rough 
surface  terrain  of  the  crystalline  regions.  Consistent  with  the  reported  improvement  of  GFA  by 
the  addition  of  Ti,  the  featureless  regions  in  the  ternary  specimen  appeared  to  be  significantly 
more  prevalent  than  in  the  graded  binary  samples. 

2.3.2  Identifying  compositional  regions  of  high  GFA  in  the  Cu-Zr-Ti  system 

The  strategy  of  applying  higher  laser  powers  stepwise  to  narrow  the  range  of  glass  for¬ 
mation  was  applied  to  deduce  the  best  glass  formers  in  the  graded  ternary  specimen.  However, 
contrary  to  the  Cu-Zr  specimens  that  featured  a  one-dimensional  gradation,  the  compositional 
gradient  in  the  Cu-Zr-Ti  specimen  was  two-dimensional,  with  every  position  on  the  specimen 
representing  a  different  composition.  Therefore,  our  experimental  method  of  obtaining  the  best 
glass  formers  required  the  fabrication  of  multiple  specimens,  each  processed  with  a  different 
laser  power  and  hence  heat  input  during  the  final  re-melting  step. 
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were  produced  according 
to  the  described  deposi¬ 
tion  scheme.  The  pro¬ 
cessing  parameters  cor¬ 
responding  to  each  spec¬ 
imen  are  specified  in  Ta¬ 
ble  3. 

Table  3  Processing  pa¬ 
rameters  for  the  laser- 
processed  Cu-Zr-Ti  com¬ 
position  libraries. 


Library  Powder  delivery  Laser  travel  speed  Final  laser 
rate(g/min)  (mm/s)  power  (W) 


Cu 

Zr 

Ti 

1 

4.12 

2.83 

0.93 

8.47  -  29.63 

200 

2 

4.12 

2.83 

1.3 

8.47  -  29.63 

240 

3 

4.12 

2.83 

1.51 

8.47  -  29.63 

280 

Measurement  of  the  complete  compositional  profile  via  EDS  for  the  three  libraries  was 
deemed  unnecessary  since  we  had  already  demonstrated  that  the  original  specimen  possessed 
the  desired  composition  gradients.  Instead,  the  approximate  boundary  positions  of  the  topo¬ 
graphically  featureless  segments  were  identified  by  DIC  imaging,  and  the  compositions  corre¬ 
sponding  to  the  endpoints  were  measured  by  EDS  to  determine  the  approximate  locations  of 
amorphous  regions  in  the  Cu-Zr-Ti  composition  map. 

The  compositional  extent  of  glass  formation  in  each  of  the  three  libraries  was  defined  by 
plotting  and  connecting  the  endpoint  compositions  to  form  enclosed  regions  on  the  Cu-Zr-Ti  ter¬ 
nary  map.  For  statistical  precision,  each  endpoint  composition  on  the  plot  represents  the  aver¬ 
age  of  six  measurements  at  the  same  position.  Figure  15  shows  the  glass  forming  region  of  all 
three  specimens  on  a  master  plot.  The  area  outlined  by  a  dashed  line  is  the  overall  composition 
space  that  was  explored.  Similar  to  the  trend  observed  in  the  graded  binary  samples,  it  is  evi¬ 
dent  that  the  amorphous  region  diminishes  with  increasing  laser  power  or  laser  line  energy. 
Based  upon  the  trend,  the  best  glass-former  within  the  explored  composition  space  was  esti¬ 
mated  to  be  Cu^jZt^jTu^.Q,  represented  by  the  black  dot  in  the  figure.  In  comparison,  the  best 
reported  glass  former  reported  by  Wang  et  al.  is  Cu57.6Zr32.4Tho,  with  a  reported  casting  thick¬ 
ness  of  4  mm.13  Inoue  et  al.  also  explored  glass  formation  in  the  Cu-Zr-Ti  ternary  system  by 
casting  and  melt-spinning  a  series  of  Cu6oZr4o-xTix  (x  =  0  -  40)  alloys.14  Cu6oZr3oTi10  was  deter¬ 
mined  as  the  best  glass  former  in  the  series,  also  with  a  maximum  casting  thickness  of  4  mm. 
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Figure  15.  A  plot  showing  the  GFA  of  Cu-Zr-Ti  alloys  within  the  Several  factors 

compositional  region  enclosed  by  the  dashed  line.  Cu51.7Zr36.7TiH.6  may  explain  the  discrep- 

was  deduced  to  be  the  optimum  glass  former.  For  comparison,  ancy  between  the  pre- 

the  red  and  blue  single  dot  indicate  compositions  reported  by  dieted  and  reported  best 

Wang  et  a/.13  and  Inoue  et  al.u  as  the  best  glass  formers  in  the  glass  formers.  First  of  all, 

ternary  system.  unlike  the  graded  binary 

Cu-Zr  specimens,  the  semi-amorphous  transition  regions  of  a  few  re-melted  lines  spanned  a 
range  of  clearly  more  than  a  few  atomic  percent,  presenting  a  greater  challenge  in  resolving  the 
endpoint  compositions  of  the  featureless  segments.  Greater  operator  errors  in  composition 
measurement  would  likely  adversely  affect  the  accuracy  of  prediction.  It  is  also  very  possible 
that  the  reported  compositions,  Cu57.6Zr32.4Tiio  and  Cu60Zr3oTi10,  are  not  actually  the  best  glass 
formers.  Inoue  and  collaborators  had  casted  their  alloys  in  large  stepwise  increments  of  5  at. 
%;  it  is  plausible  that  better  glass  formers  in  the  series  were  entirely  missed,  given  the  acute 
sensitivity  of  glass  forming  ability  to  changes  in  composition.  The  compositional  region,  with  a 
fixed  copper  content,  was  also  too  narrow  to  draw  definite  conclusions  regarding  GFA  in  the 
system.  Wang  et  a/.’s  results  were  more  convincing,  considering  that  multiple  series  of  Cu-Zr-Ti 
alloys  along  constant  Cux-Zr10o-x  tie  lines  were  interrogated.13  However,  the  ternary  tie  lines 
were  not  evenly  spaced,  with  large  compositional  gaps  between  them,  raising  doubts  regarding 
the  validity  of  their  results. 

To  validate  our  method  for  optimizing  GFA,  discrete  wedge  samples  of  Cu51.7Zr36.7TiH.6, 
Cu6oZr3oTi10,  and  (Cuo.64Zr0.36)9oTiio  were  prepared  by  copper  mold  suction  casting.  Care  was 
taken  to  maintain  consistent  processing  conditions  for  all  three  alloys.  The  wedge  critical  cast¬ 
ing  thickness  measured  by  optical  microscopy  was  1.5,  1.2,  and  1.35  mm  respectively  for 
Cu5i.7Zr3e.7Tin. 6,  Cu6oZr3oTiio,  and  (Cu0.64Zro.36)9oTiio,  confirming  the  expected  ranking  of  GFA 
among  the  three  alloys. 

Figure  16  provides  differential  scanning  calorimetry  (DSC)  heating  traces  of  the  three 
suction  cast  samples.  The  glass  transition  temperature  ( Tg ),  crystallization  temperature  (Tx), 
and  liquidus  temperature  ( 7/ )  for  each  alloy  are  listed  in  Table  2.  Also  listed  are  Trg  =  Tg/Ti  and 
AT  =  (Tx-  Tg),  the  two  most  prevalent  temperature-based  indicators  of  GFA.  It  is  noteworthy 
that  neither  parameter  accurately  predicts  the  ranking  of  GFA  among  the  three  alloys.  Accord¬ 
ing  to  AT,  Cu5i.7Zr36.7Tin. 6  is  the  poorest  glass  former.  On  the  other  hand,  the  similar  Trg  among 
the  alloys  incorrectly  suggests  that  no  differences  exist  in  their  glass  forming  ability.  These  re¬ 
sults  bring  to  reminder  the  unreliability  of  existing  temperature-based  predictors.  While  charac¬ 
teristic  temperatures  may  reflect  certain  crucial  aspects  of  glass  formation,  they  could  not  possi¬ 
bly  embody  all  the  thermodynamic  and  kinetic  factors  that  determine  an  alloy’s  ability  to  vitrify. 
Table  4  Characteristic  temperatures  and  temperature-based 
predictive  parameters  of  three  Cu-Zr-Ti  metallic  glass  alloys. 


Alloy  Composition 

Tg(  K) 

Tx(K) 

Ti(  K) 

Trg 

AT(  K) 

Cu6oZr30Tiio 

708 

751 

1146 

0.62 

43 

(Cu0.64Zr0.36)9oTiio 

693 

740 

1167 

0.59 

47 

Cu5i.7Zr36.7Ti1;L6 

682 

715 

1146 

0.60 

33 

2.4.  High  throughput  mechanical  properties  assessment  and  optimization  by  nanoinden¬ 
tation  and  dynamic  modulus  mapping 
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Although  boasting  superior  mechanical  properties  such  as  tensile  yield  strengths  ex¬ 
ceeding  1  GPa  and  the  ability  to  undergo  large  elastic  deflections  compared  to  conventional  al¬ 
loys,  monolithic  BMGs  generally  suffer  from  a  lack  of  macroscopic  ductility  in  unconfined  loading 
geometries  such  as  uniaxial  tension,  preventing  their  application  as  structural  materials.24,25 

While  crystalline  alloys  derive  their  ductility  and  toughness  from  the  abundance  and  mo¬ 
bility  of  dislocation  defects,  BMGs  on  the  other  hand,  due  to  the  complete  absence  of  disloca¬ 
tions,  cannot  accommodate  homogenous  plastic  deformation.  Instead,  plastic  strain  has  been 
shown  to  accumulate  in  the  form  of  localized  shear  bands  that  are  approximately  aligned  with 
the  plane  of  maximum  resolved  shear  stress.26,27  When  placed  under  uniaxial  tension  and 
loaded  to  failure,  monolithic  BMGs  typically  exhibit  the  growth  of  a  single  dominant  shear  band 
towards  apparent  brittle  fracture.24,28,29,30  However,  despite  their  general  reputation  as  brittle 
materials  compared  to  work-hardenable  engineering  alloys,  intrinsically  “tough”  Pd-based  and 
Zr-Cu-AI-based  BMGs  with  surprisingly  high  resistance  to  fracture  have  recently  been  discov¬ 
ered,27,31,32  prompting  the  search  for  other  monolithic  BMGs  with  high  toughness. 

The  origins  of  toughness  in  monolithic  BMG  alloys  are  not  yet  fully  understood,  but  sev¬ 
eral  empirical  criteria  based  on  elastic  constants  have  been  established.  More  specifically, 
BMG  toughness  has  been  directly  related  to  a  high  Poisson’s  ratio  (v)  or  equivalently,  a  low  ratio 
of  shear  modulus  to  bulk  modulus  (p/B). 33,34  By  plotting  the  fracture  energy  versus  v  and  p/B  for 
a  series  of  metallic  glass  alloys,  Lewandowski  et  al.  showed  that  intrinsically  tough  BMGs  must 
possess  v>  0.31-0.32  and  p/B  <  0.41-0.43  (Figure  12).33 

Elastic  constants  alone,  however,  could  not  provide  a  universal  benchmark  for  compar¬ 
ing  toughness  between  alloys.  For  example,  while  it  is  likely  that  an  alloy  with  v  =  0.4  has  high¬ 
er  fracture  toughness  than  another  alloy  with  v  =  0.2,  two  alloys  with  identical  v  are  not  neces¬ 
sarily  equally  tough.  This  is  not  surprising,  as  resistance  to  fracture  is  ultimately  dictated  by  the 
mechanisms  of  plastic,  not  elastic,  deformation  that  occur  ahead  of  the  crack  tip.  To  illustrate, 
in  2011  researchers  at  the  California  Institute  of  Technology  reported  a  damage-tolerant  Pd- 
based  BMG  that  displayed  high  fracture  toughness  values  previously  unobserved  in  metallic 
glasses.27  The  plastic  zone  size  in  fracture  toughness  tests  was  estimated  to  have  a  radius  as 
large  as  6  mm  and  contained  a  dense  population  of  shear  bands.  The  multiple  shear  band  ac¬ 
tivity  provided  a  mechanism  for  energy  dissipation  ahead  of  the  crack  tip,  leading  to  significant 
plastic  shielding  and  exceptional  frac¬ 
ture  toughness. 

More  recently,  Wang  et  al.  ex¬ 
plored  the  effect  of  surface  mechanical 
attrition  treatment  (SMAT)  on  the  ten¬ 
sile  ductility  of  BMG  tensile  speci¬ 
mens.35  Similar  to  shot  peening,  the 
SMAT  process  involves  the  repeated 
collision  of  steel  balls  against  a  speci¬ 
men  surface.  The  result  is  a  rough¬ 
ened  surface  with  a  substantial  accu¬ 
mulation  of  residual  compressive 
stress.  When  applied  to  a  Zr-based 
BMG  for  various  durations,  profuse 
shear  banding  was  observed  at  the 
surface  of  the  specimens,  with  longer 
treatments  producing  a  greater  degree 
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Figure  16.  DSC  scans  of  the  three  suction-cast 
Zr-Ti  specimens.  The  scans  were  performed  at  a 
heating  rate  of  20  K/min. 


of  shear  band  activity.  Subsequent  tensile  loading  to  failure  of  the  specimens  showed  an  in¬ 
creased  macroscopic  ductility  in  the  loading  curves  with  longer  treatment  times,  indeed  up  to  4 
%  plastic  strain.  In  contrast,  the  untreated  specimen  exhibited  nearly  no  additional  plasticity  up¬ 
on  yielding. 

These  notable  examples,  among  others,  emphasize  that  from  a  mechanistic  point  of 
view,  high  shear  band  activity  upon  yielding  of  a  BMG  alloy  is  favorable  in  terms  of  toughness 
and  ductility.  In  recent  years,  several  experimental  and  computational  studies  have  suggested 
a  direct  link  between  the  population  of  shear  transformation  zones  (STZs)  responsible  for  the 
nucleation  of  shear  bands  and  the  intrinsic  atomic-scale  elastic  heterogeneities  in  monolithic 
glasses.36,37  However,  like  GFA,  no  known  criteria  have  been  suggested  that  could  reliably  pre¬ 
dict  the  shear  band  behavior  solely  based  on  the  chemical  composition  of  an  alloy.  To  date, 
assessment  of  the  toughness  and  ductility  of  BMG  alloys  still  rely  on  conventional  ASTM  test 
methods,  which,  while  dependable,  are  incompatible  with  our  combinatorial  design  approach. 

2.4. 1_  Establishing,  compositional  trends  in  hardness  and_  modulus 

Alternatively,  instrumented  nanoindentation  provides  a  means  to  conveniently  evaluate 
the  local  (nanoscale-microscale)  mechanical  behavior  of  a  material  without  the  need  to  fabricate 
bulk  test  samples.  Unlike  more  traditional  indentation  methods,  instrumented  nanoindentation 
measures  the  load-displacement  ( P-h )  response  of  a  material,  from  which  hardness  and  re¬ 
duced  modulus  could  be  extracted.  The  hardness  is  calculated  from  the  max  load  ( Pmax )  and 
contact  depth  ( hc )  and  the  reduced  modulus  is  obtained  from  the  slope  of  the  unloading  portion 
of  the  curve.38  Deng  and  collaborators  used  instrumented  nanoindentation  as  a  combinatorial 
tool  to  explore  the  compositional  dependence  of  hardness  and  elastic  modulus  on  a  sputtered 
Zr-Cu-AI  metallic  glass  thin  film  library.39  We  have  adopted  a  similar  strategy  for  our  laser- 
deposited  graded  specimens. 

In  contrast  to  the  GFA  study,  a  high-throughput  investigation  of  the  mechanical  proper¬ 
ties  by  nanoindentation  would  benefit  from  composition  libraries  that  feature  distinct,  isolated 
compositions.  Although  lattice  mismatch  resulting  in  localized  stresses  and  strains  are  usually 
discussed  in  the  context  of  structurally  disparate  crystalline  materials,  similar  effects  could  be 
present  in  the  continuously-graded  amorphous  specimens.  As  the  glass  microstructure  varies 
across  the  compositional  profile  of  the  specimens,  the  resulting  atomic-level  strains  and  corre¬ 
sponding  stresses  could  be  manifested  as  artifacts  in  the  measured  reduced  modulus  and 
hardness  values. 

To  achieve  compositional  isolation,  Cu-Zr-Ti  compositional  libraries  were  fabricated  as  a 
square,  12  x  12  array  of  discrete  dots,  each  dot  featuring  a  fixed  chemistry  (Figure  17).  The 
dots  were  spaced  2  mm  apart.  In  the  first  step  of  the  deposition  process,  Cu  and  Zr  powder 
were  continuously  fed,  and  the  stage  followed  a  horizontal  raster  toolpath  from  the  top  left  cor¬ 
ner  of  the  array  to  the  bottom  right  corner.  To  deposit  the  hemispherical  dots,  a  250  W  laser 
was  pulsed  at  each  of  the  144  positions  for  120  milliseconds.  The  balance  of  copper  and  zirco¬ 
nium  powders  was  adjusted  from  row  to  row  such  that  a  Cu-Zr  composition  gradient  was 
achieved  along  each  vertical  column  in  the  array  while  remaining  approximately  fixed  across  a 
horizontal  row.  After  the  initial  deposit,  each  dot  was  re-pulsed  with  a  200  W  laser  for  200  milli¬ 
seconds.  Ti  was  then  deposited  with  a  150  W  laser  over  the  entire  area  of  the  initial  deposit, 
with  the  travel  speed  increasing  from  right  to  left  to  create  a  graded  Ti  profile  in  the  horizontal 
direction.  The  Cu-Zr-Ti  dots  were  re-pulsed  with  a  200  W  laser  to  homogenize  the  compositions 
and  pulsed  a  final  time  with  a  180  W  laser  for  a  short  duration  of  60  milliseconds  to  promote 
glass  formation. 
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and  hardness  values.  Each 


point  on  the  maps  represents 
the  averaged  value  of  the  four 
separate  indents.  The  fitted 
contours  displayed  smooth 
compositional  trends,  providing 
a  pallet  of  mechanical  proper¬ 
ties  whereby  alloys  could  be 
selected  for  specific  engineer¬ 
ing  applications.  The  composi¬ 
tional  trends  in  hardness  and 


modulus  also  provide  a  means 
to  further  understand  structure- 


property  correlations.  In  gen¬ 
eral,  the  mechanical  properties 
of  metallic  glass  alloys  are 
functions  of  their  chemical 


Initial  deposit 


Second  deposit 


Similar  to  the  continuously 
graded  libraries,  the  vitrified  dots  fea¬ 
tured  a  liquid-like,  smooth  surface 
that  could  be  easily  distinguished 
from  the  rough  topography  of  the 
crystallized  dots  (Figure  18).  Of  the 
144  dots,  92  were  identified  as 
amorphous.  Scanning  probe  micros¬ 
copy  of  the  vitrified  dots  yielded  sub¬ 
nanometer  surface  roughness  values 
that  were  substantially  superior  to 
mechanically  polished  surfaces.  Fur¬ 
thermore,  because  the  apex  of  a 
hemisphere  is  flat,  instrumented 
nanoindentation  could  be  performed 
directly  at  the  center  of  each  dot 
without  requiring  further  specimen 
preparation. 

To  establish  compositional 
trends  in  the  mechanical  properties 
of  Cu-Zr-Ti  glasses,  instrumented 
nano-indentation  was  carried  out 
across  each  of  the  92  compositions 
in  the  array.  Four  quasi-static  in¬ 
dents  were  performed  at  the  center  of  each  dot  in  a  fully  automated  test  sequence.  All  of  the 
high-throughput  mechanical  testing  was  accomplished  with  a  Hysitron™  Tl  950  nanoindenter, 
equipped  with  a  diamond  Berkovich  tip. 

Figure  19  shows  prop¬ 
erty-composition  maps  of  the 
measured  reduced  modulus 
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Figure  17.  (a,  top)  Photograph  of  the  Cu-Zr-Ti  dis¬ 
crete  dot  composition  library,  (b,  bottom)  Schematic 
illustrating  the  deposition  process.  Arrows  indicate  di¬ 
rection  of  increase. 


Figure  18.  DIC  optical  micrographs  showing  the  various 
surface  topographies  in  the  as-deposited  ternary  dot  libr^ 
ies.  (a-b),  amorphous  (c)  semi-amorphous,  and  (d)  crystal¬ 
line. 


bonding  character  and  the  topological  arrangement  of  the  atomic  species.  Therefore,  composi¬ 
tions  corresponding  to  local  maxima  in  the  property  maps  potentially  indicate  outstanding  char¬ 
acteristics  of  the  structure,  for  example  a  higher  degree  of  dense  atomic  packing  or  a  stiffer 
elastic  backbone. 


Hardness  (GPa)  Reduced  Modulus  (GPa) 


Figure  19.  Hardness  and  reduced  modulus  values  plotted  on  ternary  Cu-Zr-Ti  composition 
maps.  The  values  were  collected  via  instrumented  nanoindentation  measurements  performed 
on  the  discrete  dot  composition  libraries. 

2.4.2  Dynamic  rnodulus  mapping  glmonolithic  BMG  alloys 

In  recent  years,  a  number  of  computational  and  experimental  studies  have  convincingly 
demonstrated  the  spatially  heterogeneous  structure  inherent  to  amorphous  materials  and  the 
implications  of  structural  heterogeneity  on  the  mechanical  behavior  of  MGs.36'37'40'44  Egami  et 
al.  reported  that  as  much  as  25  %  of  the  volume  in  MG  alloys  is  composed  of  liquid-like  anelas- 
tic  configurations.40  Similarly,  a  notable  molecular  dynamics  study  on  a  model  Cu64Zr36  alloy 
identified  geometrically  unfavorable  soft  spots  that  exhibited  large,  non-affine  atomic-scale  de¬ 
formations.37  It  has  been  proposed  that  these  “soft”  features  are  directly  associated  with  shear 
transformation  zones  and  thus  are  responsible  for  the  plastic  behavior  of  MG  alloys.  Although 
structural  heterogeneity  has  not  been  directly  observed  experimentally,  its  existence  is  implied 
in  experimental  results  showing  the  spatial  fluctuation  of  elastic  properties.  Samwer  et  al.  used 
an  atomic  force  microscopy  indentation  method  to  probe  the  elastic  modulus  of  an  amorphous 
Pd-Cu-Si  ribbon  over  a  200  x  200  nm  area  and  compared  the  histogram  of  the  measured  values 
to  a  fully  devitrified  sample.36  The  large  Gaussian  spread  of  the  data  collected  from  the  glassy 
samples  corroborated  the  existence  of  structural  heterogeneity  in  real  MG  alloys,  and  it  was 
suggested  that  a  locally  denser  population  of  non-affine  soft  spots  are  responsible  for  the  low 
end  of  the  measured  modulus  distribution. 
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We  have  used  dynamic  modulus  mapping  in  the  nanoindenter  as  an  alternative  method 
to  probe  the  spatial  distribution  of  elastic  properties  across  the  surface  of  the  Cu-Zr-Ti  dot  spec¬ 
imens.  In  dynamic  modulus  mapping,  a  sinusoidal  load  function  with  load  amplitude  P  is  ap¬ 
plied  to  the  material  while  a  piezoelectric  actuator  scans  the  indenter  tip  over  a  user-defined  ar¬ 
ea.  The  dynamic  displacement  d  of  the  material  is  measured  as  well  as  a  phase  shift  0.  The 
storage  and  loss  stiffness  are  then  defined  as 


storage  stiffness  =  —  cos  , 

p 

loss  stiffness  =  —  sin  . 

d 

The  stiffness  values  are  easily  converted  to  modulus  values  by  applying  Hertzian  contact 
mechanics  for  a  spherical  indenter  tip.  In  all  our  mapping  experiments,  the  frequency  of  oscilla¬ 
tion  was  set  to  200  Hz  and  the  dynamic  load  amplitudes  were  chosen  to  result  in  a  1-2  nm  dis¬ 
placement  amplitude.  The  maps  were  collected  over  a  3  x  3  pm  scan  area  on  five  of  the  amor¬ 
phous  dots. 

Figure  20  provides  a  compiled  montage  of  the  storage  modulus  maps.  For  comparison 
of  the  statistical  spread  between  different  maps,  the  modulus  values  are  represented  as  a  % 
deviation  from  the  average  modulus.  Although  the  loss  modulus  maps  (not  shown)  did  not  ex¬ 
hibit  noticeable  features,  the  storage  modulus  maps  clearly  contained  distinct  regions  of  higher 
modulus  values.  From  the  well-defined  boundaries  of  the  stiff  features  in  some  of  the  maps,  it  is 
suspected  that  nanocrystals  undetectable  by  our  high-throughput  screening  method  may  be 
embedded  within  an  amorphous  matrix.  Electron  diffraction  samples  extracted  from  the  center 
of  the  dots  are  currently  being  prepared  to  confirm  whether  the  dots  are  fully  amorphous. 
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Figure  20.  Modulus  maps  collected  from  selected  dots  in  the  discrete  dot  Cu-Zr-Ti  composition 
library.  For  comparison  purposes,  the  values  are  represented  as  a  %  deviation  from  the  norm. 

3.  Explore  possible  compositions  and  composition  gradients  through  computations 

In  order  to  develop  a  full  design  methodology  for  graded  BMG  alloys  with  the  desired  gradients 
in  mechanical  properties,  we  have  developed  an  atomistic-modeling  based  design  capability  in 
parallel  to  the  experimental  efforts  in  order  to  be  able  to  form  a  validated  design  loop,  profiting 
both  from  computational  speed  as  well  as  from  experimental  validation. 

2.1  Summary:  Computational  design  capability  for  BMGs 

In  order  to  engineer  kinetic  impactors  from  graded  metallic  glass,  we  must  first  identify  sys¬ 
tems  capable  of  forming  glass,  identify  compositions  able  to  be  graded  using  LENS  or  similar 
technology,  and  estimate  the  relevant  mechanical  properties  of  the  glass.  Despite  successful 
development  of  a  few  well-understood  industrial  alloys,  identification  of  glass  forming  ability  a 
priori  is  largely  understood  through  a  set  of  empirical  rules.  Using  molecular  dynamics  simula¬ 
tions,45"47  we  have  identified  a  simple  and  computationally  efficient  method  to  describe  glass 
forming  ability  across  an  alloy  system  based  on  the  newly  proposed  “step-height  criterion”.  This 
method  has  been  tested  for  the  Cu-Zr  and  Al-La  systems  with  promising  results. 

Furthermore,  we  have  demonstrated,  as  well  known  from  previous  work,  that  MD  is  also  ca¬ 
pable  of  determining  relevant  mechanical  and  thermal  properties  making  MD  a  powerful  method 
for  identifying  and  optimizing  metallic  glass  properties. 

Finally,  in  order  to  be  able  to  model  alloys  with  many  components  and  arbitrary  constituents, 
we  have  developed  a  novel  potential  fitting  approach  that  allows  fitting  a  sensible  potential  for 
an  alloy  system  within  a  few  days.  The  new  method  is  called  Rapid  Alloy  Method  for  Producing 
Accurate  General  Empirical  Potentials”  (RAMPAGE)  and  is  freely  accessible  on  our  webpage 
atomistics.osu.edu. 

2.2  Computational  determination  of  glass  forming  ability 

In  order  to  examine  with  molecular  dynamics  (MD)  simulations  how  easily  an  alloy  can  form 
an  amorphous  glass,  several  approaches  have  been  proposed  in  the  past.  Using  a  Voronoi48 
analysis,  the  local  neighborhood  of  each  atom  can  be  tracked  as  the  system  transitions  from 
liquid  to  glass.  Atoms  in  a  liquid,49  and  subsequently  in  a  glass,  form  small  clusters  of  atoms  dic¬ 
tated  by  topological  close-packing  principles;50  clusters  cannot  efficiently  fill  space,51  promote 
amorphous  structures  and  frustrate  crystallization.  Significant  consideration  has  been  given  to 
each  atom’s  immediate  surroundings,  known  as  short-range  ordering.52  Relative  changes  in  this 
local  environment  can  be  described  in  terms  of  liquid  viscosity.  The  temperature  dependence  of 
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viscosity  is  called  “fragility”  and  has  been  used  to  investigate  glass  forming  ability  in  metallic, 
organic  and  oxide  glasses.53"55 

As  discussed  below,  packing  analyses  can  sometimes  have  misleading  results,  as  can  be 
the  case  for  fragility  modeling.  Within  this  project,  we  found  that  an  alternative  criterion,  which  is 
in  addition  computationally  rather  efficient,  can  give  more  reliable  predictions  of  GFA.  The  data 
presented  in  Figure  21  shows  the  fraction  of  atoms  with  a  coordination  number  of  12  as  a  func¬ 
tion  of  temperature.  Coordination  12  was  chosen  because  it  directly  competes  with  the  FCC  and 
HCP  crystal  structures.  Large,  steep  changes  are  poor  glass  formers,  while  small  steps  tend  to 
be  good  glass  formers  as  shown  in  Figure  21.  As  shown  in  Figure  22,  the  step  height  may  cor¬ 
relate  with  liquid  fragility  in  Al-La.  We  peopose  that  this  correlation  could  hold  the  key  to  a  priori, 
computation-based  determination  of  glass  forming  ability. 

Determining  fragility  requires  a  separate,  expensive  simulation  that  can  represent  half  or 
more  of  the  total  computational  expense  for  a  given  alloy  composition.  Even  after  expending  this 
effort,  the  resulting  data  can  be  unreliable,  requiring  statistical  analysis  to  use  in  any  scientifical¬ 
ly  rigorous  way.  Finding  the  step  height  requires  data  already  stored  as  part  of  the  liquid-to- 
glass  quench  simulation.  If  step  height  could  be  definitively  proven  to  correlate  to  fragility,  it 

could  reduce  the 
computation  time 
required  for  each 
alloy,  while  eliminat¬ 
ing  the  least  reliable 
data  provided  by 
simulations.  The  ex¬ 
perimental  determi¬ 
nation  of  fragility  re¬ 
quires  sensitive 
measurements  on  a 
supercooled  liquid56 
or  extrapolation  of 
high  temperature 
data  to  the  glass 
transition.  By 
contrast,  an  experi¬ 
ment  to  find  the  step 
height  would  require 
measurements  of 
the  liquid  structure 
at  high  temperature 
and  the  final  glass: 
measurements  that 
have  already  been 
made  consistently.57 
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Cu65Zr35  is  a  glass  former.  Note  the  gradual 
change  prior  to,  and  in  the  supercooled  regime. 


Cu50Zr50  is  a  bulk  glass  former.  Note  the  near-zero 
structural  change  from  liquid  to  glass. 
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AI70La30  is  not  a  glass  former.  Note  the  sudden 
transition  deep  in  the  supercooled  regime. 


AI70La30  is  a  glass  former.  Note  the  gradual 
transition  deep  in  the  supercooled  regime. 


Figure  21:  Fraction  of  atoms  with  coordination  number  12  as  a  function 
of  temperature  shown  for  four  glass  forming  simulations:  (a)  Cu65Zr35, 
(b)  Cu5oZr5o,  (c)  AI7oLa3o,  (d)  AI30La7o.  Good  glass  formers  (a),  (b)  and 
(d),  have  small  step  heights  and  gradual  transitions  from  liquid  to  glass. 
Poor  glass  formers,  (c),  have  tall,  sudden  steps. 
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Without  experimental  input,  this  method  can  qualitatively  rank  candidate  compositions,  and 
possibly  bound  the  glass  forming  range.  With  the  inclusion  of  some  experimental  data  from  Ref. 
58,  separation  of  compositions  capable  of  forming  glass  becomes  possible.  As  shown  in  Figures 
22  and  24,  no  Al-La  alloy  forms  glass  when  the  step  height  exceeds  0.105.  Using  this  threshold, 
the  gap  in  glass  forming  from  La13  to  La43  can  be  successfully  filtered.  This  finding  is  also  con¬ 
sistent  with  the  understanding  that  high  fragility  can  indicate  poor  glass  forming  ability.  Such  a 
threshold  may  exist  in  Cu-Zr,  but  the  narrow  band  of  fragilities  displayed  in  that  system  does  not 
provide  sufficient  resolution  to  verify  it  (Fig.  23). 
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Figure  22:  Relationship  between  liquid  fragility  and  step  height 
for  the  Al-La  glass  system.  The  rough  correlation  between  the 
two  could  form  the  basis  for  reliable  computation-first  determi¬ 
nation  of  glass  forming  ability,  as  fragility  is  the  best  available 
indicator  of  experimental  glass  forming  ability.  If  this  relation 
ship  holds  for  other  systems,  it  could  eliminate  both  a  computa¬ 
tional  end  experimental  trouble-maker. 


Figure  23:  Step  heights  for  all  simulated  Cu-Zr  compositions.  Ex¬ 
perimental  glass  formers  are  green  and  experimental  non-glass 
formers  are  black.  The  marker  indicates  the  fraction  of  atoms  with 
coordination  number  12  in  the  liquid,  and  the  arrow  points  to  the 
concentration  of  coordination  number  12  in  the  glass.  Cu5oZr5o 
had  no  discernable  step  height  and  is  also  an  experimental  bulk 
glass  former,  capable  of  forming  amorphous  sections  over  1  cm 
thick.  No  threshold  in  step  height  filters  the  glass  forming  section 
due  to  crystallization  in  the  simulation  above  Zr70. 


AILa  Liquid-Glass  Step  Summary 
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Figure  24:  Step  heights  for  all  simulated  Al-La  compositions.  Ex¬ 
perimental  glass  formers  are  blue  and  experimental  non-glass 
formers  are  black.  The  marker  indicates  the  fraction  of  atoms  with 
coordination  number  12  in  the  liquid,  and  the  arrow  points  to  the 
concentration  of  coordination  number  12  in  the  glass.  Al-La  is 
unique  in  that  two  discrete  glass  forming  ranges  exist  with  a  gap 
spanning  La17  to  La43.  This  gap  can  be  successfully  filtered  as 
shown  in  Figure  22.  The  existence  of  the  gap  would  likely  preclude 
grading  for  compositions  near  its  ends. 

2.3  RAMPAGE  interatomic  potentials 

The  step-height  method  proposed  in  Sec.  2.2  shows  promise  for  binary  alloys,  but  the 
most  deployable  metallic  glasses  have  3  or  more  components.  In  order  to  simulate  these  com¬ 
plex  systems  with  MD  requires  matching  interatomic  potentials.  Multi-component  potentials  can 
be  produced,  but  always  trade  accuracy  for  ease  of  production.  The  easiest  method,  the  John¬ 
son  Alloy  Model  (JAM),59  linearly  interpolates  the  elemental  potentials  according  to  the  electron 
density  term.  This  method  is  exceptionally  efficient  but  has  been  shown  to  cause  unphysical  re¬ 
sults,  particularly  in  amorphous  systems.  The  traditional  method  of  creating  highly  accurate  po¬ 
tentials  takes  a  holistic  view  of  the  system  and  aims  to  incorporate  as  many  test  structures  as 
possible.60  A  single  potential  can  take  months  or  years  to  generate  with  this  approach. 

In  order  to  test  many  metallic  glass  alloys,  a  tool  was  created  to  produce  potentials  su¬ 
perior  to  the  JAM,  but  with  substantial  automation  over  the  traditional  method.  The  “Rapid  Alloy 
Method  for  Producing  Accurate  General  Empirical  Potentials”  (RAMPAGE)  was  developed  to 
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quickly  fit  multi-component  Finnis-Sinclair-type  EAM  potentials  from  DFT  and  experimental  data. 
Using  this  tool,  the  glass-forming  method  previously  described  could  be  extended  to  actual  en¬ 
gineering  applications  in  metallic  glass.  Creating  and  supporting  this  tool  has  required 

As  an  impactor  penetrates  a  medium,  it  encounters  significant  friction  and  substantial 
heating.  Critical  metallic  glass  properties,  like  viscosity  and  density,  change  rapidly  as  tempera¬ 
tures  approach  and  exceed  the  glass  transition  temperature  (Tg).  Increases  in  temperature  can 
lead  to  crystallization,  melting  or  complete  destruction  of  the  glass.  Therefore,  thermal  proper¬ 
ties  and  the  expected  change  in  mechanical  properties  should  be  considered  for  impactor  de¬ 
sign.  MD  simulations  track  changes  in  viscosity  of  the  glass  with  temperature  and  the  glass 
transition  temperature  can  be  found  with  reasonable  experimental  agreement  as  shown  in  Fig¬ 
ures  25  and  27  for  CuZr,  and  Figures  26  and  28  for  AILa.  MD  employs  a  very  fast  quench  rate, 
slightly  elevating  the  simulated  glass  transition  temperature.61,62  Coupled  with  FEM  calculations, 
this  information  can  determine  the  ideal  dimensions  to  ensure  a  payload  will  survive  to  a  particu¬ 
lar  depth. 
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Figure  25:  Simulated  and  experimental  glass  transition  tempera¬ 
tures  for  Cu-Zr  metallic  glass  alloys.  Experimental  data  from  Refs. 
63  and  64. 


Al-La  Glass  Transition  Temperatures 
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Figure  26:  Simulated  and  experimental  glass  transition  temperatures 
for  Al-La  metallic  glass  alloys.  Experimental  data  from  Ref.  66. 
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Cu50Zr50  Volume  on  Cooling 
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Figure  27:  Volume  per  atom  as  a  function  of  temperature.  The  change  in 
volume  with  temperature  (dV/dT)  is  fit  to  a  polynomial.  The  derivative  of  the 
polynomial  maximizes  at  the  glass  transition  temperature.  As  an  impactor 
heats  due  to  friction,  changes  in  volume  per  atom  will  affect  density  and 
strength  which  were  found  to  be  primary  design  criteria.  Simulations  can  de¬ 
termine  the  maximum  operating  temperature  for  a  given  alloy  and  provide 
the  dependence  of  density  on  temperature. 
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X  MD  Simulation  ■  Experiment  [Li  et  al.  2008]  ■  Experiment  [Calvayrac  et  al.  1983] 

Figure  28:  Experimental  and  simulated  mass  density  for  Cu-Zr.  Exper¬ 
imental  data  from  Refs.  67  and  68. 

2.4  Computational  prediction  of  metallic-glass  properties 

Molecular  dynamics  provides  other  necessary  design  information.  As  shown  in  Figure 
28,  MD  reproduces  the  mass  density  of  Cu-Zr  metallic  glass  very  well.  As  the  FEM  results  indi¬ 
cate  that  density  is  a  powerful  design  parameter  this  makes  MD  a  powerful  tool  for  identifying 
metallic  glasses  best  suited  to  impactor  applications. 

Molecular  dynamics  can  provide  qualitative  estimates  of  a  material’s  yield  strength  from 
simulated  stress-strain  curves  as  shown  in  Figure  29.  Due  to  the  high  strain  rates  simulated  in 
MD,  the  simulated  yield  strength  and  Tg  are  overestimated  but  compositional  trends  and  sys¬ 
tem-wide  behaviors  are  preserved.69,70  The  same  simulation  that  provides  yield  data  also  pro¬ 
vides  the  moduli.  Coupled  with  FEM  predictions  discussed  in  Sec.  1,  these  additional  design 
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Cu-Zr  Amorphous  Density 
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parameters  extracted  from  MD  simulations  can  be  used  to  optimize  impactor  materials  design, 
validated  by  the  experiments  described  in  Sec.  2. 


Simulated  Cu-Zr  Yield  Strength 
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Figure  29:  Simulated  yield  strength  of  Cu-Zr  glasses.  The  simulated  strain 
rate  is  lO^/sec.  Due  to  this  high  strain  rate,  the  yield  strengths  presented 
here  are  overestimates.  However,  the  general  compositional  trend  should  be 
preserved.71,72 
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Appendix  -  Simulation  Details  and  Constitutive  Description  of  Materials 

1 . 1  Units  used  in  Abaqus 


Length 

Meter 

Force 

Newton 

Time 

Second 

Mass 

Kg 

Density 

Kg/m3 

Stress 

N/m2 

Young’s  Moduli 

N/m2 

1 .2  Part  and  Assembly 

1.2.1  Metallic  Projectile 
See  figure  4 

1.2.2  Sand 


Material:  Sand 
(Eulerian  Part) 

Section:  10x10 
Length:  50 

The  Eulerian  part  can  be  divided  into  several 
layers  with  void,  full  or  discrete  material. 


1.2.3  Concrete  Bunker 


Material:  Concrete  (Lagrangian  Part) 

In  order  to  make  the  sand  part  are  endless,  there 
is  a  concrete  confinement  outside  the  lateral  side 
of  sand  part. 

Outer  section:  12x12 

Inner  section:  10x10 

Length:  5 1 

Thickness:  1 


1 . 3  Material  Properties 

1.3.1  Concrete 

1.3.1. 1  Elastic,  Thermal  and  Electrical  Properties 

The  material  properties  used  for  concrete  bunker  are  come  from  examples  in  ABAQUS  example 
manual.  Properties  of  concrete  are  same  as  Koyna  dam,  which  is  a  classic  example  used  for 
seismic  analysis. 


Density 

2643 

Elastic  Moduli 

31.02E+09 

Poisson  ratio 

0.2 

Conductivity 

0.7 

Specific  Heat 

880 

1.3. 1.2  Concrete  Damage  Plasticity 

(1)  Plasticity 


Dilation  Angle 

Eccentricity 

fto/fco 

K 

Viscosity  parameter 

36.31 

0.1 

1.15 

0.0007 

O(default) 

(2)  Compressive  Behavior 


Yield  Stress 

Inelastic  Strain 

1.30E+07 

0 

2.41E+07 

0.001 

(3)  Tensile  behavior  (Type:  Displacement) 


Yield  stress 

Displacement 

Damage  Parameter 

2.90E+06 

0 

0 

1.94E+06 

6.62E-05 

0.381217 

1.30E+06 

0.00012286 

0.617107 

8.73E+05 

0.000173427 

0.763072 

5.86E+05 

0.00022019 

0.853393 

3.92E+05 

0.000264718 

0.909282 

2.63E+05 

0.000308088 

0.943865 

1.76E+05 

0.00035105 

0.965265 

1.18E+05 

0.000394138 

0.978506 

7.92E+04 

0.000437744 

0.9867 

5.31E+04 

0.000482165 

0.99177 

1.3.2  Sand 

1.3.2.1  Elastic,  Thermal  and  Electrical  Properties 

Most  of  the  general  properties  of  sand  are  estimated  from  the  average  value  of  dry  sand. 


Density 

5000 

Shear  Moduli 

5.0E+07 

Conductivity 

0.2 

Specific  Heat 

830 

1. 3.2.2  Equation  of  State 

The  equation  of  state,  compaction  state,  are  especially  coming  from  Brown,  J.  L.  et  al.  (2007). 
Linear  us—up  relationship  is  used  for  sand,  and  the  Gruneisen  parameter,  T0,  is  assumed  to  be 
1.0  in  this  model. 


us—up  Liner  Relationship 

tis  < 

:0  +  sup  =  0.243  +  2.348up 

Reference  speed 

Slope 

Gruneisen  parameter 

(c0,  m/s) 

(s) 

(To) 

243 

2.348 

1.0 

P  —  a  EOS  Compaction 

Reference  sound  speed  in  the  porous  material ,  ce,  m/s 

145.8 

Value  of  the  porosity  of  the  unloaded  material,  n0 

0.4 

Pressure  required  to  initialized  plastic  behavior,  pe,  Pa 

0.9E+06 

Compaction  pressure  at  which  all  pores  are  crushed  ,  ps.  Pa 

1.13E+09 

1.3.3  Metallic  Glass  (Zr55Al10Ni5Cu30) 

Most  of  the  material  properties  of  metallic  glass,  Zr55Al10Ni5Cu30,  are  determined  by  Mashimo, 
T.  et  al.  (2006).  Some  undetermined  parameters  are  instead  of  common  glass  properties  because 
of  the  absence  of  experiments.  The  plasticity,  failure  model,  and  the  equation  of  state  are  also 
assumed  from  experimental  data. 

1.3.3. 1  Mechanical  Properties 

1.3.3.2  Elastic,  Thermal  and  Electrical  Properties 


Density 

6790 

Elastic  Moduli 

85.7E+09 

Poisson  ratio 

0.37 

Shear  Moduli 

31.28E+09 

Conductivity 

5.00 

Specific  Heat 

700 

1.3.3.3  Johnson  -  Cook  Plasticity 

Johnson-Cook  Plasticity  is  suitable  for  high-strain  rate  deformation  of  most  metals,  and  most 
applicable  to  dynamic  failures.  Johnson-  Cook  plasticity  model  belongs  to  Mises  plasticity 
model,  and  must  be  in  conjunction  with  either  the  linear  elastic  material  model  or  equation  of 
state.  In  this  simulation,  the  Mie-Gruneisen  equation  of  state  material  model  is  selected  to 
establish  the  hydrodynamic  material  model. 

The  stress  of  metallic  glass  can  be  assumed  to  be  the  following  equation. 

a  =  [A  +  B  (e^)n]  [1  +  C  In  (e)]  [1  +  T m] 

where  sPl  is  the  equivalent  plastic  strain  and  A,  B,  n  and  m  are  material  parameters  measured  at 
or  below  the  transition  temperature.  C,  is  the  material  parameter  depend  on  strain  rate,  and  i  is 
strain  rate.  The  temperature  dependence  term  is  defined  as  the  following  equation. 

(  0  forT  <  Ttr 

f  =  \(T-  Ttr)/(Tm  -  Ttr)  ,  for  Ttr  <  T  <  Tm 
l  1  for  T  >  Ttr 

where  T  is  the  non-dimensional  temperature,  T  is  the  current  temperature,  Tm  is  the  melting 
temperature,  and  Ttr  is  the  transition  temperature. 

In  Abaqus  explicit,  Johnson-Cook  plasticity  can  be  conjunction  with  the  progressive  damage  and 
failure  model,  which  based  on  equivalent  plastic  strain.  Johnson-Cook  damage  model  can  also 


depend  strain  rate  and  temperature  as  well.  The  dependencies  are  assumed  to  be  the  forma  of  the 
following  equation. 

£fl  =  [di  +  d2Exp(d3n)][l  +  d4lni][  1  +  d5f] 


Where,  dl5  d2,  d3,  d4,  and  d5  are  failure  parameters.  In  Abaqus  simulation,  the  elements  that 
meet  the  failure  criterion  are  deleted. 

The  material  parameter  used  in  this  simulation  is  displayed  in  the  following  table. 


Johnson-Cook  Plasticity  Model 

A 

B 

n 

C 

m 

324E+06 

2890E+06 

0.56 

0.0025 

1.09 

Johnson-Cook  Damage  Model 

d i 

d2 

d3 

d 4 

ds 

0.54 

4.78 

4.47 

0.014 

1.12 

2.5 

2 

1.5 

v: 

v: 

0 

LA 

i 

0.5 

0 

( 

Johnson-Cook  Plasticity 

* 

* 

I_ I_ I_ I_ I_ l_ I 

)  0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4 

Strain 

1.3.3.4  Gruneisen  Equation  of  State 

The  tabular  function  values,  A  and  f2,  are  calculated  from  results  from  Hugoniot  compressive 
test  in  Mashimo,  T.  et  al.  (2006). 

Pfavol)  Afevoz)  T  A 

The  data  of  Hugoniot  compressive  test  are  shown  in  table  below.  The  Gruneisen  parameter,  T0,  is 
assumed  to  be  1.5  which  is  identical  to  common  glass  temperately. 


Hugoniot  compressive  test 

Specimen 

thickness 

Initial 

density 

Impact 

velocity 

Shock 

velocity 

Particle 

Velocity 

Pressure 

Density 

No. 

(m) 

(m/s) 

(m/s) 

(m/s) 

(m/s) 

(GPa) 

(Kg/m3) 

1 

1.816E-03 

6790 

670 

5260 

171 

6.11 

7020 

2 

2.082  E-03 

6750 

758 

5300 

185 

6.62 

6990 

3 

1.810  E-03 

6760 

1033 

5360 

160 

5.8 

6970 

4 

1.819  E-03 

6750 

1169 

5290 

380 

13.57 

7270 

5 

1.825  E-03 

6740 

1193 

5230 

390 

13.75 

7280 

6 

1.831  E-03 

6760 

1081 

5139 

767 

26.65 

7950 

7 

1.826  E-03 

6750 

1547 

5255 

1103 

39.12 

8540 

8 

1.828  E-03 

6740 

1756 

5400 

1248 

45.42 

8770 

1.3.3.5  Drucker  -  Prager  Plasticity 

Shear  criterion:  exponent  form  (developed  from  Mohr-coulomb  yield  surface) 

(1)  Dilation  angle 

Dilation  angle,  ip,  measured  in  the  p-q  plane  at  high  confining  pressure,  unit  in  degree. 

The  angle  of  dilation  controls  an  amount  of  plastic  volumetric  strain  developed  during 
plastic  shearing  and  is  assumed  constant  during  plastic  yielding.  The  value  of  \|/  =  0 
corresponds  to  the  volume  preserving  deformation  while  in  shear. 

(2)  Flow  potential  eccentricity 

Flow  potential  eccentricity,  e,  defines  the  eccentricity  between  yielding  function  and  its 
asymptote. 

(3)  Yield  surface 

The  Drucker-Prager  plasticity  is  established  on  yield  surface. 


According  to  the  confining  stress,^,  and  loading  stress,  cr3,  cq  =  a2  =£  cr3,  the  yield 
surface  of  material  can  be  established  by  finding  a  best  fit  curve.  Equation  2.1  is  the 
exponential  function  of  yield  surface  for  a  material. 

F  =  aqb  —  p  —  pt  =  0  ( o.i ) 


The  material  parameters  a,  and  b  can  be  given  directly,  or  subscribed  from  triaxial  test 
alternatively,  a  =  —  b  =  In  this  simulation,  the  value  of  a  and  b  are  estimated 

from  the  Mohr  Coulomb  yield  surface  of  Zr55Al10Ni5Cu30. 


1.3.3.6  Damage:  Ductile  Damage  for  Metal 


The  triaxial  test  (or  compressive  and  tensile  tests),  the  flyer  plate  test,  and  ballistic  penetration 
test  are  required  to  build  the  JH-2  damage  model  for  a  material.  The  Hugoniot  state,  pressure, 
and  some  common  material  properties  of  metallic  Glass  ( Zr55Al10Ni5Cu30 )  can  be  determined 
by  Mashimo  T.  et  al.  (2006).  Holmquist,  T.J.  et  al.  (1995)  described  the  method  and  results  of 
establishing  the  JH-2  model  of  float  glass.  Combine  with  Drucker-Prager  Model  and  Mie 
Gruneisen  EOS,  the  JH-2  model  can  calibrate  in  ductile  damage  model  in  Abaqus.  JH-2  model  is 
suitable  for  brittle  material,  progressive  damage,  large  pressure  and  high  strain  rate.  All  the 
parameters  used  in  JH-2  model  can  be  determined  from  Mashimo  T.  et  al.  (2006)  and  Lund,  A.  C. 
and  Schuh,  C.  A.  (2003).  Equations  2.2  and  2.3  show  the  mathematical  relationship  between 
constant  pressure,  P,  fracture  strain,  s1’1,  stress  triaxiality,  q,  and  strain  rate,  s.  The  continuous 
relationship  between  fracture  strain,  stress  triaxiality,  and  strain  rate  can  be  determined  according 


to  the  Hugoniot  test  from  Mashimo  T.  et  al.  (2006),  which  offers  eight  results  of  flyer  plate  tests 
with  similar  strain  rate,  which  can  determine  a  continuous  relationship  between  stress  triaxiality, 
strain  rate,  and  fracture  strain. 
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'  intact 
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&HEL 
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PiO 


+  ■ 
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[1  +  C7n(£(i))] 


(0.3) 


yP HEL  '  PhEL / 

where,  Phel  is  pressure  at  Hugoniot  elastic  limit;  is  stress  at  Hugoniot  elastic  limit;  T  =  pt 
Three  viable  are  required  in  ABAQUS  shown  below. 


Fracture  Strain 

i£Pl) 

Stress  Triaxiality 

in) 

strain  rate 

is?1) 

Round  off 

i£pl) 

0.36579 

-21.2643 

88398 

1.00E+05 

0.32387 

-19.7921 

88857 

1.00E+05 

0.23781 

-16.4196 

94163 

1.00E+05 

0.14233 

-11.7558 

208906 

2.00E+05 

0.14093 

-11.6758 

213699 

2.00E+05 

0.08453 

-7.89775 

418897 

5.00E+05 

0.08016 

-7.53763 

604053 

5.00E+05 

0.07748 

-7.31234 

682713 

5.00E+05 

Note:  strain  rate  inputted  in  ABAQUS  are  all  assumed  to  be  identical  for  enough  information. 


Metallic  Glass  iZr3BAl10Ni5Cu30) 

Density,  p0 

6750Kg/m3 

Young’s  Moduli,  E 

85.7GPa 

Poisson  Ratio,  v 

0.37 

Shear  Moduli,  G 

31.28GPa 

Strength  Model 

A 

0.6027 

N 

0.5393 

C  (consider  with  different  stain  rate) 

- 

Normalized  intact  strength:  of  =  AiP*  +  T*}w(l  +  C  In  e) 

Tensile  strength,  at 

1.53Gpa 

Compressive  strength,  oc 

1.77Gpa 

HEL 

11.07  GPa 

Pressure  at  HEL,  PHEL 

P  =  Kl[i  +  K2p2+K3p3 

4.87  GPa 

Stress  at  HEL,  aHEL 

6.2GPa 

Volume  strain  at  HEL,  l*iHEL 

,=f-i 

Po 

0.0392 

Drucker-Prager  Plasticity 

3 

Yield  Surface:  2.5 

F  =  aqb  —  p  —  pt  =  0, 

1.5 

1 

V  =  —  ^  (2ai  +  03). 

0.5 

q  =  (ai  ~  cr3)  0 

* 

L 

■* 

.  /X 

ti¬ 

lt 

2-10  1  : 

p 

> 

Phel 

(Phel^)1^ 

0.4223 

1 

b  =  — 

N 

1.8452 

.  1 

pt  =  a<rcD  +-crc 

1.8073  Gpa 

Dilation  angle,  ip 

O(default) 

Flow  eccentricity ,  e 

0.1  (default) 

Equation  of  State  :  P(.£voi)  =  fi(.£Voi )  +  /2  Pm 

(Bulk  Moduli) 

110.8  GPa 

480.6  GPa 

K3 

-1201  GPa 

Pressure,  P  ( GPa ) 


Damage  Model  (use  model  of  glass) 

Dx  (glass)  0.043 

D2  (glass)  0.85 


1.4  Mesh  (Explicit) 


(1)  Element  type  used  for  metallic  projectile  is  C3D6.  (Lagrangian  solid  part) 
Distortion  control=0.1;  Element  deletion=yes;  Max  Degradation^. 5  (max  is  1); 
Seed  size=0.025  Element  shape=Wedge; 

(2)  Element  type  used  for  concrete  bunker  is  C3D8.  (Lagrangian  solid  part) 
Distortion  control=0.1;  Element  deletion=yes;  Max  Degradation^ ; 

Seed  size=0.25;  Element  shape=Hex; 


(3)  Element  type  used  for  sand  is  C3D6  (Hex).  (Eulerian  part) 

^  Element  Type  |  S3  | 


Element  Library 

Family 

Standard  ®  Explicit 

Geometric  Order 

o  Linear  Quadratic 

Hex  [wedge  [  Tet 


0  Thermally  coupled 
Element  Controls 

Second-order  accuracy;  ©  Yes  (®)  No 

Hourglass  control:  (®>  Use  default  ©  Relax  stiffness  ©  Stiffness  ©  Viscous  ©  Combined 

Stiffness-viscous  weight  factor  0.5 

Scaling  factors;  Displacement  hourglass:  |l^  Linear  bulk  viscosity;  1  Quadratic  bulk  viscosity:  1 


EC3D8R:  An  8-node  linear  eulerian  brick,  reduced  integration,  hourglass  control. 


1 .5  Predefined  Fields  and  Boundary  Condition 
1.5.1  Initial  Condition 


The  projectile  has  an  initial  velocity  of  600m/s  in  z-direction.  The  Eulerian  part  is  consisted  of 
void  parts  and  sand  parts. 


r  Edit  Predefined  Field  1^1 

Name:  Predefined  Field-sand 

Type:  Material  assignment 

Step:  Initial 

Instance:  base-1 
Volume  Fractions 


Definition:  l@>  jUniformj  0  Discrete  fields 


Region 

base-1. 

sand-1 

Void 

1 

base-1.  sand_l 

0 

1 

2 

base-1. sand_2 

0 

1 

3 

base-1. sand_3 

1 

0 

4 

base-1. sand_4 

1 

0 

Note:  When  the  same  region  is  specified  more  than 
once,  the  last  entry  is  used. 


OK  |  Cancel  | 


Initial  velocity  of  projectile 


Predefined  field  of  Eulerian  part 


1.5.2  Boundary  Condition 


The  surrounding  of  concrete  model  is  confined  in  three  directions  translation  and  rotation. 


Cancel 


Edit  Boundary  Condition 

sal 

Name:  BC-1 

Type:  Displacement/Rotation 

Step:  Step-1  (Dynamic,  Explicit) 

Region:  (Picked) 

CSYS:  (Global)  ; 

Distribution:  Uniform 

^  f(x) 

[V]  Ul:  [o 

1 

(V]  U2:  [o 

[7]  U3:  [o 

1 

[7]UR1:  0 

radians  1 

m  UR2:  [o 

radians  1 

[7]  UR3:  [o 

radians  1 

Amplitude: 


(Instantaneous) 


Note:  The  displacement  boundary  condition 
will  be  reapplied  in  subsequent  steps. 


The  inflow  surface  and  outflow  surface  of  Eulerian  part  are  applied  with  Eulerian  Boundary. 


•v  Edit  Boundary  Condition 
Name;  BC-3 
Type:  Eulerian  boundary 

Step:  Step-1  (Dynamic,  Explicit) 

Region:  (Picked)  [;> 


r  Edit  Boundary  Condition 
Name:  BC-5 
Type:  Eulerian  boundary 

Step:  Step-1  (Dynamic,  Explicit) 

Region:  (Picked)  [ 


Flow  type:  o  Inflow  Outflow  (_)  Both 


Flow  type:  Inflow  o  Outflow  Both 


Inflow:  Free 


Outflow:  Equilibrium 


Cancel 


Cancel 


Inflow  Surface 


Outflow  Surface 


1.6  Additional  Results 


1.6.1  Radially  Graded  Properties 


The  following  figures  show  effects  of  bullet  with  graded  properties  along  thickness.  The  remark 
“1”  in  the  figure  denotes  to  larger  properties,  and  “3”  denotes  to  smaller  properties.  The  depth  of 
each  situation  is  listed  in  the  following  table. 


Original 

Density  - 

Density+ 

Stiffness- 

Stiffness+ 

Strength- 

Strength+ 

36.5196 

35.9109 

36.9396 

36.5799 

36.8868 

36.8891 

37.8805 

- 

-1.64% 

+1.15% 

+0.17% 

+1.00% 

+1.01% 

+3.73% 

velocity  velocity 


displacement 


Graded  Properties  along  thickness 


0  5  10  15  20  25  30  35  40 


displacement 


1.6.2  Concrete  Penetration 


Simulations  were  performed  with  concrete  slabs  at  a  depths  of  20,  25,  30  and  40m.  The 
homogenous  projectile  was  able  to  penetrate  the  20  and  25m  deep  slabs  but  not  the  30m  deep 
slab.  The  projectile  stopped  in  the  sand  before  reaching  the  40m  deep  projectile. 


U-V 
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